Entropically-stabilised growth of a two-dimensional random tiling 
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The assembly of molecular networks into structures such as random tilings and glasses has recently 
been demonstrated for a number of two-dimensional systems. These structures are dynamically- 
arrested on experimental timescales so the critical regime in their formation is that of initial growth. 
Here we identify a transition from energetic to entropic stabilisation in the nucleation and growth 
of a molecular rhombus tiling. Calculations based on a lattice gas model show that clustering 
of topological defects and the formation of faceted boundaries followed by a slow relaxation to 
equilibrium occurs under conditions of energetic stabilisation. We also identify an entropically- 
stabilised regime in which the system grows directly into an equilibrium configuration without the 
need for further relaxation. Our results provide a methodology for identifying equilibrium and non- 
equilibrium randomness in the growth of molecular tilings, and we demonstrate that equilibrium 
spatial statistics are compatible with exponentially slow dynamical behaviour. 



The properties of two-dimensional supramolecular net- 
works have been the focus of growing interest in re- 
cent years with most efforts directed towards the con- 
trolled introduction of translational order into such sys- 
tems [U [2] . However there have been several recent ob- 
servations of surface-bound supramolecular arrays which 
assemble into dynamically-arrested structures akin to 
glasses [3-5 which lack translational order. Such ar- 
rangements raise many interesting questions related to 
the growth of random systems [6, . In particular it is 
important to distinguish randomising effects which arise 
from kinetic effects, such as nucleation [8[|9], from equi- 
librium disorder due to entropic terms in the free energy. 
Entropically-stabilised disorder may be regarded as in- 
trinsic randomness, whereas kinetically-driven disorder 
is often determined by sample history and preparative 
conditions. In one recent study [3 a random molecu- 
lar rhombus tiling was shown to have equilibrium (max- 
imum entropy) spatial correlations, despite being frozen 
on an experimental timescale. In such a system the max- 
imum entropy configuration must form, and be frozen 
in, during the initial growth, since the spatio-temporal 
fluctuations which normally facilitate the evolution of 
kinetically-trapped configurations to equilibrium are ab- 
sent. However, it is not clear, a priori^ that there is a set 
of local rules for molecular attachment which can lead to 
the direct growth of a 'perfect', i.e. maximum entropy, 
configuration. 

In this paper we address this question and show that 
equilibrium and non-equilibrium effects in the growth of 
a rhombus tiling [1017 may be distinguished using tile- 
tile correlations of arrays simulated using a lattice gas 
model [18]. Direct growth to a configuration with equi- 
librium statistics occurs when entropic terms dominate 
the free energy, while non-equilibrium effects result in 
faceted islands and clustering of topological defects. 

The parameters which control growth are the tile-tile 
interaction energy, £, the tile adsorption energy, /i, and 
the temperature, T. We consider (see Fig. 1) a trian- 
gular lattice with sites, labelled that are either occu- 




10"^ 10° 10^ 10^ 10^ 10' 

t [MCS] 



0.00** — ' — ' — ' — ' — ^ 



1.7 



1.9 



2.0 



FIG. 1: Schematic of: (a) lattice and tiles; (b) defect diffusion 
mediated by tile detachment and re- attachment; (c) annihila- 
tion of a defect pair (left to right) or generation of defect pair 
(right to left), (d) The fraction of empty sites as a function 
of time (in units of Monte Carlo sweeps) (ksT = 0.3, lattice 



size N = 10 . (e) Dependence of Cg 
text and in (d)] on /x. 



and Ac [defined in 



pied by half a tile or vacant. Each rhombus tile occupies 
two adjacent sites and lies in one of three orientations 
(distinguished by different colours). For molecules de- 
posited from solution [3] jn corresponds to the difference 
between the net adsorption and solvation energies of the 
molecule (/i > implies a preference for solvation). The 
binding energy per tile for a completely tiled surface is 

Ebind = -(2£ - /i). 

Using a Metropolis algorithm [19^, sites are chosen 
randomly and, if empty, a tile is added with probabil- 
ity e"^^/^^^ for AE > and with unit probability for 
AE < 0, where AE is the associated change in energy. If 
the site is occupied, tile removal is accepted with prob- 
ability e-^^/^^^/3. The factor of 1/3 ensures detailed 
balance is satisfied. The energy required to remove a tile 
with p nearest neighbours is Erem = {p^ — m), which, 
for /i > 0, may be either positive or negative depend- 
ing on the local environment. One Monte Carlo sweep 




(MCS) corresponds to the random inspection of 3A^ sites 
of the (rhomboid) lattice {N is the maximum number 
of adsorbed tiles), and sets the unit of time. Periodic 
boundary conditions are used and all energies are hence- 
forth expressed in units of e. This is a generalisation, 
through the introduction of the parameter, /i, of a model 
previously used to show that rhombus tilings are glassy 

m- 

As the time increases, the fraction of empty lattice 
sites, c(t), reduces from 1 (empty lattice) and eventually 
relaxes to an equilibrium phase at a /i-dependent con- 
stant value of c(t), Ceq. In the initial growth phase c{t) 
falls until an abrupt change in gradient occurs. At this 
point, to a good approximation, there are no more avail- 
able vacancies (neighbouring pairs of unfilled triangles) 
which could directly accommodate a tile. However, the 
lattice is not completely tiled and triangular void defects 
are also present. These are topological defects with two 
effective charges corresponding to triangles pointing up 
and pointing down [Ml [iQl [H] • Further relaxation is me- 
diated by defect diffusion and annihilation (neighbouring 
defects of opposite effective charge form a vacancy which 
may be occupied by a tile - see Fig. 1) . Since defect 
diffusion is an activated process (barrier 3 — /i) there is a 
slowing down which gives rise to the clear change in gra- 
dient discussed above. In the equilibrium regime there is 
a dynamic balance between the generation of triangular 
defect pairs (from the removal of tiles) and their diffusion 
and annihilation ^18^ . 

We parameterise the value of c{t) after the initial 
growth phase as Cg (see Fig. 1 where the values for Cg 
and Ceq are identified for the /i = 1.7 curve). We also 
introduce a parameter, Ac(= Cg — Ceq)^ to quantify the 
difference between the defect density in equilibrium and 
immediately after the initial growth phase. The depen- 
dence of the parameters c^, Cgg, and Ac on ji is shown 
in Fig. 1(e) over the parameter range 1.7</i<2.1. As 
expected, Ceq increases with increasing /i since the energy 
barrier for tile detachment is reduced. Interestingly, in 
the range ji > 2 the binding energy, Ei)ind > 0, and no 
tiling would be expected for an ordered system. However, 
random tilings do grow in this regime; the variation of 
c{t) for /i = 2 is shown in Fig. 1(d), and values for c^, 
Ceq and Ac extracted in the regime where £^5ind > 
(~ 2.1 > fi > 2), are shown in Fig. 1(e). We show 
below that in this regime entropic contributions lead to 
a free energy, F, given by {NEi^ind — TS) where S is 
the entropy, which can be negative, favouring a tiling, 
even when Eijind > (Joseph et al [22 make similar 
arguments in the context of entropically-stabilised qua- 
sicrystals [131 US]). Furthermore, in this entropically- 
stabilised regime, Ac — > 0, implying, as confirmed below, 
that the initial growth phase leads directly to an equilib- 
rium regime. 

We now consider the differences in nucleation, mor- 
phology and tile statistics in the energetically-stabilised 




FIG. 2: (a-c) Simulated growth of tilings for varying /x with 
ksT = 0.2: (a)/i = 0.5, high nucleation density, irregular 
islands (c(t) = 0.7); (b) /i = 1.5, reduced nucleation den- 
sity, facetted islands (c(t) = 0.7); (c) /x = 1.8, example of 
an inhomogeneous tile distribution within strongly facetted 
island; (d) schematic of growth along a straight interface for 
1 < /i < 2. 

(/i < 2) and entropically-stabilised (/i > 2) regimes, fo- 
cussing first on low values of /i. In Fig. 2 we show islands 
which have been nucleated and are growing in the initial 
growth regime. For < /i < 1, Eremi the barrier for tile 
removal, is positive even for p = 1 indicating that any nu- 
cleated island formed by two neighbouring tiles is stable 
(note for /i < even isolated tiles are stable nuclei). Ac- 
cordingly, the simulated growth in this regime [Fig. 2(a)] 
shows a large number of small, irregular islands of tiles. 
This is an essentially homogeneous growth regime: very 
quickly the islands merge forming an imperfect tiling of 
the plane. 

For 1 < /i < 2 we have heterogeneous growth; as ji 
increases, islands become larger and their number de- 
creases. Furthermore, the islands that form are faceted 
and hexagonal [see Fig. 2(b)] with a clear deficit of tiles 
in one of the three possible orientations (colours) in each 
of the six triangular segments of the island [see Fig. 2(c)]. 
The faceting occurs since the smallest energetically stable 
nucleus requires a minimum of three tiles in a hexagonal 
configuration. Outward growth results in the (approxi- 
mate) propagation of the hexagonal shape since growth 
along an edge favours the addition of a row formed by one 
of the two tile orientations with an edge parallel to the 
island boundary. For example, in Fig. 2(d), rows of blue 
and/or red tiles grow and where they meet an upward- 
or downward-pointing triangular defect is formed. The 
downward defect is trapped, but the upward defect may 
be occupied by a green tile, leading to an excess of de- 
fects of one effective charge (downward-pointing in this 
case) in each segment. Inspection of the hexagonal is- 
land in Fig. 2(c) shows rows of tiles of two colours with 
either defects, or a tile of the third colour where they 
meet, consistent with this simple explanation. Thus the 
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faceting is due to minimisation of the boundary energy, 
but also results in clustering of defects with the same 
effective charge. A local imbalance of tiles gives rise to 
an increase in entropic free energy, and is not expected 
for an equilibrium configuration. We note the interesting 
similarity between the growth of energetically-favoured 
hexagonal tilings observed here with the 'arctic circle' 
problem in rhombus tilings subject to hexagonal confin- 
ing boundaries [151 ES] • 

As growth in this regime continues the islands merge 
while maintaining, approximately, the primordial struc- 
ture introduced in the nucleation stage. An example is 
shown in Fig. 3(a) which shows the tiling that is formed 
immediately after all the growing islands merge. Defects 
in this tiling are not distributed uniformly as in equilib- 
rium [17]. Fig. 3(b) shows the topological charge density 
for the tiling of Fig. 3(a). Defect clustering is evident 
in the large variations in topological charge. The origi- 
nal nucleation sites within the tiling can be identified as 
singular points in this defect density [Fig. 3(b)], which 
confirms the spatial correlation of defect clustering and 
nucleation sites. Further temporal evolution governed 
by defect diffusion and annihilation [18 leads to equili- 
bration of the tiling, see Fig. 3(c), when c{t) reaches its 
equilibrium value [see Fig. 3(d)]. 

The tilings are analysed using a lifting dimension 
[T3l [24] in which an effective height, h{f) is assigned 
to each vertex (with in-plane co-ordinates, r) in the 
tiling. The height is calculated using the scheme shown 
in the Fig. 3(e) inset in which a displacement along a 
rhombus edge leads to a change in height of ±1. The 
height correlation function, C{r) = {[h{0) — /i(r)]^), 
can be calculated and for a maximally random tiling, 
C{r) = (7ri^o)~^ ln(r) + c, has a logarithmic dependence 
on position, where c is a constant and Kq = 7t/9 [13^ . 
Fig. 3(e) shows the correlation functions during the sim- 
ulated growth of the tilings in Figs. 3(a) & 3(c). For 
the tiling in Fig. 3(a) the correlation function is not 
logarithmic. For increasing times the correlation func- 
tions approach a linear dependence on ln(r), with the 
expected gradient 9/7r^, confirming that the final con- 
figuration [Fig. 3(c)] is equilibrated. This supports the 
hypothesis that a logarithmic dependence is associated 
with an equilibrium configuration rather than kinetically- 
controlled randomness, but the exponentially slow ap- 
proach to equilibrium cannot account for tilings which 
are both dynamically-arrested and maximum-entropy. 

For /i > 2 nucleated islands do not show faceting or in- 
homogeneities [Fig. 4(a) inset] and our simulations show 
direct growth into a maximum-entropy configuration. In 
Fig. 4(a) we plot correlation functions for tilings immedi- 
ately after the initial growth regime is completed [deter- 
mined by the change in gradient in c{t)] and find an ap- 
proach to a logarithmic dependence on r as /i increases. 
Note from Fig. 1(e) that Ac ^ for /i > 2 and these 
results confirm that this simple parameter provides a re- 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

In (r) 

FIG. 3: (a) An inhomogeneous tiling resulting from strongly 
faceted growth of multiple islands, ksT = 0.2, /a = 1.8, 
AT = 1.6 X 10^ (- 4 X 10^ tiles shown), (b) Clustering of 
topological charge shown in a charge density map. The value 
at a point corresponds to the number of upward- pointing mi- 
nus downward-pointing defects within a range of three times 
the average defect separation, (c) The tiling after relaxing 
to an entropically-maximized equilibrium state (^ 4 x 10^ 
tiles shown), (d) the corresponding c(t) behaviour, and (e) 
height correlation functions during relaxation, indicating con- 
vergence to maximum randomness. 



liable indicator for a regime of direct growth into an equi- 
librium configuration without the requirement for defect- 
mediated relaxation. 

To confirm that tilings for /i > 2 are entropically sta- 
bilised we need to establish the equilibrium phase bound- 
ary /i(T) between the tiled and empty phases. We estab- 
lish this by investigating whether an interface between 
an equilibrated tiling and an empty lattice recedes (no 
growth) or propagates (growth) The value of /i 

where this transition occurs is plotted against temper- 
ature in Fig. 4(b). As discussed above, for an ordered 
system no growth is expected for /i > 2, but the free en- 
ergy may be negative when Ebind > 0, if Suie > Eund/T, 
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FIG. 4: (a) Height correlation functions calculated for tilings 
immediately after the initial growth stage [where c(t) = Cg] 
at ksT = 0.3 with varying /i. (b) Calculations of stationary 
interfaces marking the tiled-untiled phase equilibrium, (c) 
Illustrative representation of the /a — T parameter space, in- 
dicating the region where tilings are entropically- stabilized. 



where Stue is the entropy per tile, or /i < 2-\-TStiie' Our 
simulations [Fig. 4(b)] give Stue = in excellent 

agreement with the ideal value for the rhombus tiling 
entropy density of 0.323066A:s [131125]. 

These results are highly relevant to recent experiments 
[3]. We propose that, in both the energetically- and 
entropically- stabilised regimes, growth proceeds through 
the initial regime to the point identified in Fig. 1 where 
there is a slowing down of the evolution of the tiling. 
Further evolution is determined by the barrier to defect 
propagation, (3 — /i). If this is small compared with 
/c^T, equilibration can occur through defect propaga- 
tion. However, for many molecular systems the barrier is 
at least an order of magnitude greater than the thermal 
energy and the configuration is therefore dynamically- 
arrested with spatial statistics which are frozen immedi- 
ately after the initial tiling of the surface. A broad range 
of possible spatial distributions can occur, including, for 
/i ^ 2, the recently-observed maximum-entropy arrange- 
ment [5. However for lower values of /i a configuration 
with a frozen-in, non-equilibrium spatial distribution of 
tiles [comparable to Fig. 3(a)] might be attainable in ex- 
periments. 

The phase behaviour observed in these simula- 



tions is summarised in Fig. 4(c) and invites anal- 
ogy with magnetic Ising systems since the total en- 
ergy of a partially tiled surface is given hy E = 

-f Ei5i (lfj=i ^3 - 1) + 2 Ei5i where rii = 1(0) 
for an occupied (unoccupied) site (for a related example 
of the application of the lattice gas model to adsorbed 
molecular layers see Ref. [26]). The index i runs over all 
triangular sites and j runs over the three nearest neigh- 
bours of site i. The tile-tile interaction, e, is analogous to 
the spin-spin coupling, normally denoted by J, fi is anal- 
ogous to magnetic field, and rii to spin state. In Fig. 4(c) 
a boundary at /i = 2 shows the threshold above which the 
internal energy is positive. However ioi T < Tc ^ £ /ks 
there is an entropically-stabilised regime for the random 
rhombus tiling [shaded region in Fig. 4(c)]. This phase 
boundary can be determined in our simulations up to 
ksT 0.5 which we identify as an approximate criti- 
cal temperature for this transition. The deposition of a 
molecular layer is thus equivalent to a quench from high 
ji (analogue magnetic field). The growth dynamics af- 
ter such a quench allows an investigation of this phase 
diagram for systems where dynamics are slow. 

Our results show that equilibrium and non-equilibrium 
randomisation may be distinguished for the rhombus 
tiling. Moreover we have shown that equilibrium spatial 
statistics may occur even for dynamically-arrested sys- 
tems although other outcomes such as faceting and defect 
clustering are also possible. These results have general 
relevance for molecular layers adsorbed at a liquid-solid 
interface where it has previously been assumed, correctly 
in many cases, that a dynamic equilibrium is established 
with molecules continually exchanged between solvated 
and adsorbed states (see [27^ for example). Our results 
show that such exchange is not required for the forma- 
tion of maximum-entropy arrangements. There are also 
interesting links between entropically-stabilised growth 
and several other problems in bio- and condensed matter 
physics, such as the crystallisation of anisotropic parti- 
cles ^28j. Furthermore, molecular rhombus tiles provide 
a new system to explore, both experimentally and the- 
oretically, equilibrium and non-equilibrium behaviour in 
connection with 'Coulomb' and other exotic phases which 
can exhibit fractional excitations [24 , such as frustrated 
magnets with effective magnetic monopoles [29 , quasi- 
crystals [7 , and glasses [6| [30]. 
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